Time-Resolved Measurement of a Charge Qubit 
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We propose a scheme for monitoring coherent quantum dynamics with good time-resolution and 
low backaction, which relies on the response of the considered quantum system to high-frequency ac 
driving. An approximate analytical solution of the corresponding quantum master equation reveals 
that the phase of an outgoing signal, which can directly be measured in an experiment with lock-in 
technique, is proportional to the expectation value of a particular system observable. This result 
is corroborated by the numerical solution of the master equation for a charge qubit realized with a 
Cooper-pair box, where we focus on monitoring coherent oscillations. 

PACS numbers: 42.50.Dv, 03.65.Yz, 03.67.Lx, 85.25. Cp 



An indispensable requirement for a quantum computer 
is the readout of its state after performing gate opera- 
tions. For that purpose it is sufficient to distinguish be- 
tween two possible logical states. At the same time, it is 
desirable to demonstrate the coherence of time evolution 
explicitly. For solid-state qubits, this has been accom- 
plished by Rabi-type experiments [H, @] . In general the 
qubit state measurement is destructive, so that an inter- 
ference pattern emerges only after a number of experi- 
mental runs. Certainly, it would be preferable to observe 
signatures of coherent dynamics already in a single run. 

Both the charge and the flux degree of freedom of su- 
perconducting qubits can be measured by coupling the 
qubit to a low-frequency "tank" circuit that is excited 
resonantly [H, 0| ■ In doing so one makes use of the fact 
that the resonance frequency of the slow oscillator de- 
pends on the qubit state which, in turn, influences the 
phase of the oscillator response [!, 0, HI . The drawback of 
this scheme, however, is that the coherent qubit dynam- 
ics is considerably faster than the driving. Thus, one can 
only observe the time-average of the qubit state, but not 
time-resolve its dynamics. Measuring the qubit by driv- 
ing it at resonance is possible as well This however 
induces Rabi oscillations, making the qubit dynamics dif- 
fer significantly from the undriven case 0, Hf ■ 

Here, by contrast, we propose to probe the qubit by a 
weak high-frequency driving that directly acts upon the 
qubit without a tank circuit being present. We find that 
the resulting outgoing signal possesses sidebands which 
are related to a phase shift and demonstrate that the 
latter is related to a qubit observable. Validating this 
relation numerically for a Cooper-pair box, we show that 
the underlying measurement scheme principally enables 
monitoring the coherent qubit dynamics experimentally 
in a single run with good time-resolution and fidelity, 
whereas the backaction on the qubit, induced by the driv- 
ing, stays at a tolerable level. 

Dissipative quantum circuit. — Although later on we fo- 
cus on the dynamics of a superconducting charge qubit 
as sketched in Fig. [IJ our measurement scheme is rather 
generic and can be applied to any open quantum system. 



We employ the system-bath Hamiltonian [^, [13, U | 
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where Hq denotes the system Hamiltonian and Q is a 
system operator. The bath is modelled by LC circuits 
with charges qk and conjugate momenta pk, where Ck 
and Lk are effective capacitances and inductances, re- 
spectively, and Afe are the corresponding coupling con- 
stants. It is convenient to introduce the spectral density 
I(uj) = | J2k Afc(Cfc/ifc) 1 / 2 (5(w — uJk) which we assume 
to be ohmic, i.e. I(uj) = wZq with an effective impedance 
Zq [ll|, [T3, EH- By standard techniques we obtain the 
Bloch-Redfield master equation for the reduced system 
density operator p in the weak-coupling limit [LH ]. 

p = ~{H ,p] - ~[Q, [Q,p]]-ij±[Q, [Q, />]+], (2) 



with the anticommutator [A, _B]_ 
operators Q = i [Hq ,Q]/h and 



AB + BA and the 



dr / duj S(lo) cos(lot)Q(— t). (3) 
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Here, S(ui) = J(w) coth(htu /21cbT) is the Fourier trans- 
form of the symmetrically ordered equilibrium correla- 
tion function ^ ([£in(i"), £in(0)]+) oq at temperature T with 
regard to the collective bath coordinate £; n = X^fcAfeQfe- 
The notation X(t) is a shorthand for the Heisenberg op- 
erator Uq (t)XUo(t), where Uq is the system propagator. 
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FIG. 1: (color online) Cooper pair box (CPB) coupled to a 
transmission line with ohmic effective impedance Zq. 
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In order to relate the quantum dynamics of the cen- 
tral circuit to the response via the transmission line, we 
employ the input-output formalism which is an es- 
tablished tool in quantum optics and has also been used 
in quantum circuit theory [5j. It starts from the Heisen- 
berg equation of motion for the environmental mode k 
which reads q k + u\q k = u)\.Q, where uj k = (L k C k )~ 1/2 
denotes the angular frequency of mode k. Owing to its 
linearity, this equation of motion can be solved formally. 
Inserting the obtained solution into the Heisenberg equa- 
tions of motion for the system operators, one arrives at 
the so-called quantum Langevin equation T(| 17 1. For an 
ohmic environment, the latter possesses the inhomogene- 
ity £m{t) — Z Q(t), where the noise operator £ in is fully 
determined by the correlation functions given above. 

Alternatively, one can write the quantum Langevin 
equation in terms of the outgoing fluctuations 15]. The 



result differs only by the sign of the dissipative term, so 
that the inhomogeneity now reads £ ou t(i) + Z Q{t). The 
difference between these two equations relates the input 
and the output fluctuations via 



£out — £in — —2ZqQ — — 
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[H ,Q], 
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which is a cornerstone of the input-output formalism [151 ] 
and holds for any weakly coupled £j n . We used that for 
weak dissipation, Q rss i\Ho,Q]/H is essentially bath- 
independent. 

Response to high-frequency driving. — We next probe 
the system by driving it via the transmission line with an 
ac signal Acos(ilt) that also couples to the system op- 
erator Q. Then, the Hamiltonian acquires an additional 
term: Ho — > Ho + QAcos(flt), and the master equa- 
tion Q changes accordingly. For the input £j n , this cor- 
responds to one coherently excited incoming mode such 
that (£in(i)) = Acos(ttt), while the r.h.s. of the input- 
output relation Q remains unchanged. 

Because the driving must not significantly alter the 
system dynamics, we assume that the amplitude A is 
sufficiently small, so that the driving can be treated per- 
turbatively. This yields the ansatz p(t) = po(t) + pi(t), 
where po(t) is the unperturbed state. To lowest order 
in A, pi obeys pi = CoPi — j{A{Q,po} cos(fii), where Co 
denotes the superoperator on the r.h.s. of Eq. ©. This 
linear inhomogeneous equation of motion can be solved 
formally in terms of a convolution between the propaga- 
tor of the undriven system and the inhomogeneity. If the 
driving frequency is much larger than all relevant sys- 
tem frequencies, one may separate time scales to obtain 

p(t) = p (t) - i^[Q, Po (t)} sin(m) , (5) 

which identifies eA/Ml as the necessarily small pertur- 
bation parameter. Together with the input-output rela- 
tion (0| this solution allows us to compute the response 
of the system. 



In an experiment it is possible to employ a lock-in tech- 
nique with the incoming signal providing the reference os- 
cillator. This singles out the high-frequency components 
of the outgoing signal, which correspond to the second 
term of the density operator (0 and read 

1AZ 

(Ut(t))=Acos(ilt)+-^ ([[H ,Q],Q}} ahx(Qt), (6) 

where (. . .)o = tr[po(i) . . .] refers to the undriven dynam- 
ics. Writing next (£out(*)) = A! cos [fit — <^ f (t)], we find 
the central expression 



[[H0,Q],Q]) : 



(7) 



which relates a small phase shift <$}f(t) between the in- 
put and the output signal to a hermitian system observ- 
able describing the unperturbed low-frequency system 
dynamics. This means that the time- resolved evolution 
of the open quantum system can be monitored in a sin- 
gle run by continuously measuring the phase shift <^ f (t) 
with appropriate experimental techniques. In this con- 
nection, Eq. |(7J) represents the basis for our proposed 
measurement scheme. Below we will explore its valid- 
ity and limitations for a specific system by comparing 
the phase of the output (£out(t)) with the expectation 
value ([[i?Oj Q]i Q])qi both computed from the numerical 
solution of the master equation ^ in the presence of ac 
driving. 

Monitoring coherent qubit dynamics. — A particular 
case of a quantum circuit which recently attracted much 
interest is a Cooper pair box (CPB) which is sketched in 
Fig. [1] and described by the Hamiltonian 



H, 



CPB 
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where N is the number of excess Cooper pairs in the 
box, so that the charge operator reads Q — 2eN = 
2e J2n N\N)(N\. The charging energy Ec is determined 
by various capacitances, while the scaled gate voltage N g 
and the effective Josephson energy Ej are controllable. 
If the charging energy is sufficiently large and N g ~ 1/2, 
only the two charge states |0) = ||) and |1) = ||) matter 
and form a qubit [j], 0, [ll[ described by the Hamiltonian 



H, 



qb _ 



(9) 



where the Pauli matrices Oi are defined in the qubit sub- 
space and E e i = 4i?c(l — 2Ag). The qubit energy split- 
ting reads Sw q b = (E^ + E 2 ) 1 / 2 . Moreover, Q qh = ea z 
while by virtue of relation (J7J) the phase of the output 
signal is linked to the qubit observable o~ x according to 



(*) 



h 2 n 



(10) 



3 



1 

0.5 


-0.5 
-1 




- r 



'Mi 't f» A « . 



— (e(0> • • < f (0 



won - 



10 



20 
t [h/Ej] 



30 



40 



10° 




2A£2 






2co qb 
-«-»► 




10~ 2 


-— |<e(co)>| 2 








■ --- |(e qb (co))| 2 






io- 4 


<»qb 

' 4 


£2 

4 





5 10 20 
co \Ej/h] 



50 



FIG. 2: (color online) Decaying qubit oscillations with ini- 
tial state |t) in a weakly probed CPB with 6 states for 
a = Z e 2 /h = 0.08, A = 0.1Ej/e, £ c = 5.25£j and 
iV 9 = 0.45, so that E cl = 2.1Ej and u qb = 2.3£j/ft. (a) Time 
evolution of the measured difference signal (Q) oc (£out) — (£in) 
(in units of 2eEj/h) of the full CPB and its lock-in amplified 
phase <f> ou t (frequency window Af2 = 5Ej/h), compared to 
the estimated phase {&x)o m the qubit approximation. 

The inset resolves the underlying small rapid oscillations with 
frequency £1 = 15-Bj /ft in the long-time limit, (b) Power spec- 
trum of (Q) for the full CPB Hamiltonian (solid) and for the 
two- level approximation (dashed). 



This means that the high-frequency component of (Q), 
which is manifest in the phase of the outgoing signal Q , 
contains information about the low-frequency qubit dy- 
namics in terms of the unperturbed (ct x }q. 

We now turn to the question how relation (fTT)|) allows 
one to retrieve information about the coherent qubit dy- 
namics in an experiment. Figure [21a) shows the time 
evolution of the expectation value (Q(t)) for the initial 
state ||) = |1), obtained via numerical integration of the 
master equation @ with the full Cooper-pair box Hamil- 
tonian |8]) in the presence of the ac driving which in prin- 
ciple may excite higher states. The driving, due to its 
rather small amplitude, is barely noticeable on the scale 
chosen for the main figure, but only on a refined scale for 
long times; see inset of Fig. [21a). This already insinuates 
that the backaction on the dynamics is weak. In the cor- 



responding power spectrum of (Q) depicted in Fig. [2b), 
the driving is nevertheless reflected in sideband peaks at 
the frequencies Q and f2±cj q h>. In the time domain these 
peaks correspond to a signal cos[f2i — 4> ut{t)]- More- 
over, non-qubit CPB states leads to additional peaks at 
higher frequencies, while their influence at frequencies 
lu < f2 is minor. Experimentally, the phase <j) ou t{t) can 
be retrieved by lock-in amplification of the output signal, 
which we mimic numerically in the following way [l8| : We 
only consider the spectrum of £ out in a window f2 ± Af2 
around the driving frequency and shift it by — Q. The 
inverse Fourier transformation to the time domain pro- 
vides <pout(t) which is expected to agree with <^^f(t) and, 
according to Eq. (fTU)) , to reflect the unperturbed time 
evolution of (cr x )o with respect to the qubit. Although 
the condition of high-frequency probing, Q 3> w q b, is 
not strictly fulfilled and despite the presence of higher 
charge states, the lock- in amplified phase <j> ut(t) and the 
predicted phase (p^At) are barely distinguishable for an 
appropriate choice of parameters as is shown in Fig.[^a). 

In order to quantify this agreement, we introduce the 
measurement fidelity F = (4> on t, (cx)o)) where (/,<?) = 
J dt fg/(J dt f 2 J dt g 2 ) 1 / 2 with time integration over the 
decay duration. Thus, the ideal value F = 1 is assumed 
if 4>ont(t) ancl ( cr x(i))o are proportional to each other, 
i.e. if the agreement between the measured phase and 
the unperturbed expectation value (a x )o is perfect. Fig- 
ure [3Ja) depicts the fidelity as a function of the driving 
frequency. As expected, whenever non-qubit CPB states 
are excited resonantly, we find fcl, indicating a signif- 
icant population of these states. Far-off such resonances, 
the fidelity increases with the driving frequency f2. A 
proper frequency lies in the middle between the qubit 
doublet and the next higher state. In the present case, 
il w 15Ej/h appears as a good choice. Concerning the 
driving amplitude, one has to find a compromise, because 
as A increases, so does the phase contrast of the outgo- 
ing signal (JB]) , while the driving perturbs more and more 
the low- frequency dynamics. For the frequency chosen 
above, the fidelity is best in the range A — 0.1-lfiw q b/e. 
This corresponds to eA/Kl w 10 _2 -10 -3 , which justifies 
our perturbative treatment. 

Measurement quality and backaction. — For any quan- 
tum measurement, one has to worry about backaction on 
the system in terms of decoherence. In our measurement 
scheme, decoherence plays a particular role, because both 
the driving and the ohmic environment couple to the 
CPB via the same mechanism. This is reflected by the 
fact that the predicted phase ([TO]) is proportional to the 
dimensionless dissipation strength a — Zoe 2 /h. However, 
a should not exceed 0.1 in order to preserve a predom- 
inantly coherent time evolution, which also means that 
our measurement is weak, but destructive. The condi- 
tion a < 0.1 together with the above conditions on the 
driving amplitude and frequency provides phase shifts 
4> ou t of the order 1°, which is small but still measurable 
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FIG. 3: (color online) (a) Fidelity defect 5F = 1 - F and 
(b) time-averaged trace distance between the driven and the 
undriven density operator of the CPB for various driving am- 
plitudes as a function of the driving frequency. All other 
parameters are the same as in Fig. [2] 



system coupled to a dissipative environment. The cru- 
cial requirement for this is the possibility to drive the 
system coherently via one environmental high-frequency 
mode accompanied by measuring the phase of the outgo- 
ing signal via lock-in techniques. By analyzing the high- 
frequency response, we have found that the phase of the 
output signal is related to a particular system observable. 
We have substantiated this relation by computing both 
quantities numerically for a charge qubit implemented 
with a Cooper-pair box. For decaying coherent oscilla- 
tions, we have demonstrated experimental feasibility on 
condition that the coupling to the environment is not too 
weak and that the driving frequency exceeds the qubit 
splitting and is off resonance with higher levels. Then 
the measurement fidelity is rather good, while the low- 
frequency qubit dynamics is almost not affected by the 
driving and transitions to higher levels do not play a rele- 
vant role. The implementation of our scheme will enable 
the demonstration of quantum coherence of solid-state 
qubits in single-shot experiments. 

This work has been supported by the DFG through 
SFB 631 and by the German Excellence Initiative via 
"Nanosystems Initiative Munich (NIM)" . 



with present technologies. The additional decoherence 
due to the driving, by contrast, is not noticeable. This 
is in agreement with the first-order result tr(p 2 ) = tr(pQ) 
which follows from Eq. ([5]). 

In order to investigate to what extent the driving af- 
fects the quantum state of the CPB, we compute the trace 
distance D(t) — |tr|p(<) — po(t)\ [l9[ between the den- 
sity operators of the driven system p(t) and the undriven 
reference po(t). Its time average D quantifies the per- 
turbation due to the driving. Figure EJb) indicates that 
D oc A/fl, unless the driving is in resonance with higher 
levels. This confirms the picture drawn by studying the 
measurement fidelity F. For practically all parameters 
used in Fig. [3j we found that the total population of 
levels outside the qubit doublet is always less than 0.1%. 
The only exception occurs again in the case of resonances 
with non-qubit CPB levels. Far from these resonances, 
the system is faithfully described with only the qubit lev- 
els. 

In our investigations, we have not considered excita- 
tions of quasiparticles which are relevant once the driv- 
ing frequency becomes of the order of the gap frequency 
of the superconducting material. Thus, for an aluminum 
CPB our model is valid only for f2 < 100 GHz. Since a 
typical Josephson energy is of the order of some GHz, a 
driving frequency f2 w 20Ej/h is still within this range 
while it provides already a good measurement quality. 

Conclusions. — We have proposed a method for the 
time-resolved monitoring of the dynamics of a quantum 
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